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Abstract 

This study examined one of the most used techniques in 
Structural dynamics testing - Experimental Modal Analysis 
(EMA). This analysis determines structural modal parameters 
i.e. natural frequencies, damping factors and the 
characteristic mode shapes of the predominant modes of 
vibration, from measured transfer function data. 

Once these dynamic parameters were known, they could be 
compiled into a modal model of the structure. Modifications 
in the form of additions (or removal) of mass, stiffness, or 
damping could then be implemented analytically to predict 
the dynamic behavior of the modified structure. This 
mathematical technique called Strucural Dynamics 
Modification (SDM), allows for the analytical evaluation of 
the dynamic effects of modifications thereby dispensing with 
experimentation with actual physical modifications. 

The first half of this study involved the validation 
and necessary software development for calibrated EMA and 
SDM. In the second half, EMA and SDM were combined and their 
accuracy together was evaluated. This first involved the 
construction of a modal model from a calibrated EMA on a 
test structure. Then using the modal model and SDM, 
predicted modification results were obtained and compared to 
results from an EMA on the corresponding physically modified 
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1. INTRODUCTION 


1.1 FEA,EMA and SDM 

An important consideration in the design and operation 
of mechanical structures and machines is a knowledge of 
their dynamic characteristics. These characteristics are 
determined by the interaction of the mass, elastic and 
damping properties of the system and can be expressed 
through modal parameters. These modal parameters include the 
damping, natural frequencies and mode shapes of the 
predominant modes of vibration and can be found either 
analytically or experimentally. Once known, the modal 
parameters can be compiled to form a complete system dynamic 
model or modal model. This model then allows the prediction 
of the dynamic response of the structure or machine due to 
any applied set of forces. 

Simple systems can be described by continuous models 
where the modal parameters are derived from the formulation 
of one or more differential equations. Often however, the 
complexity of the system is such that approximations such as 
discretization of a continuous structure is used. In this 
case, the actual system is modelled by a set of discrete 
masses or rigid bodies coupled by springs and dampers. 
Application of the basic equations of dynamics will then 
result in a system of second order differential equations 


which can be solved simultaneously for the modal parameters. 
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More recently for complex structures, the dynamics have 
been calculated by means of Finite Element Analysis (FEA) 
(for example Zienkiewicz [32]). FEA involves the 
discretization of a mechanical structure into small 
elements. The structures dynamic equations are then 
assembled as a combination of the individual elements 
equations of motion and boundary conditions. To accurately 
model the structure, many elements and thus hundreds or 
thousands of degrees-of-freedom may be required. Because of 
the size and complexity of these models, they can be 
expensive and time consuming to both develop and use. In 
addition, it is often found that due to an insufficient 
number of elements and unrepresentative boundary conditions 
that the model does not accurately describe the actual 
Situation. Because of these problems, validation of 
analytical models is often performed using experimental 
dynamic testing. 

One of the most popular forms of dynamic testing in 
recent years is Experimental Modal Analysis (EMA). For a 
detailed explanation of the subject see Richardson [24]. 
This analysis is based on obtaining simultaneous excitation 
and response measurements at points of interest on the 
structure. By dividing the Fourier transform of the response 
by the Fourier transform of the excitation, a transfer 
function between the various response and excitation points 
is formed. The modal parameters are then extracted by curve 


fitting an analytical expression to the experimentally 
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determined transfer functions. The modal parameters can then 
be compiled to form an experimental modal model. 

Experimental modal models have several advantages over 
the analytical. The number of degrees-of-freedom (m) in the 
finite element model dictates the models size of m x m. 
Therefore, m, or hundreds to thousands of simultaneous 
equations must be solved. Although the solution results ino 
modes, the large majority of them are outside the frequency 
range of interest. In the experimental case, the curve fit 
results in n modes in the required range. Thus the model 
Size can be considerably smaller (n x n). This means that 
for the experimental model, computation is faster, memory 
requirements are less, and thus a smaller computer can be 
used. The second advantage is that an experimental model 
represents the actual structure. The third is that an 
experimental model can be obtained for structures too 
complex to be easily modelled analytically. 

The combination of analytical and experimental 
analysis, when possible, represents a most powerful tool in 
assessing a structures dynamics. By direct comparison 
between the experimental and analytical modal model, 
shortcomings Or errors in each can be detected. It then 
becomes possible to either re-model portions of the finite 
element model or re-measure data at particular points on the 
structure: torform’ a’setvof more’ consistent results: 

Once the modal model is acceptable, the effects of 


proposed changes to the structure to either optimize the 
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dynamic behavior or alter it to alleviate problems can be 
examined usSing a technique called Structural Dynamics 
Modification (SDM). This technique developed by engineers at 
Structural MeaSurement Systems Inc., is outlined in some 
detail in Formenti and Welaratna [6]. SDM uses either an 
experimental or analytical model to analytically evaluate 
the effects of modifications to the structure in the form of 
stiffeners, dampers or point masses. The technique is very 
concise and thus numerous modification proposals can be 
investigated quickly. This fact takes on immense importance 
especially in FEA. Because of the complexity of the 
structure usually modelled, placement and amount of 
modification is not obvious and thus the evaluation of 
modifications can involve a fair amount of trial and error. 
To alter and run a finite element model for each change can 
be time consuming and expensive. However, by first finding a 
Salient change quickly and easily using SDM, the finite 
element model could then be run to further detail the effect 


of the particular modification. 


1.2 Thesis Outline 

The fundamental objective of this study was to extend 
the basic EMA software package developed by Fyfe [7] to 
enable SDM. This involves further software development as 
well as transfer function calibration and the introduction 


of the SDM technique. 
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Chapter two summarizes modal analysis and the 
relationship between theory and experiment. Also included is 
the method of extracting modal parameters from transfer 
function data. Further theoretical investigations extend the 
concepts to the SDM technique and its development using the 
modal model. 

Chapter three discusses the practical aspects of using 
EMA and SDM. These considerations include the relevant 
fundamentals of signal analysis and hardware requirements in 
relation to data acquisition and data processing. This leads 
iPromavcisccussicnoorthne curve tatting mrocessmand transter 
Function calibration. AS well, practical considerations of 
the SDM technique involve the use of the Newton-Raphson 
technique in determining the modified structure eigenvalues. 

The fourth chapter discusses the results obtained from 
the separate validation and combined evaluation of 
calibrated EMA and SDM. As the accuracy of the predicted 
dynamic effects of modifications on a structure using SDM 
depends on the accuracy of modal parameter estimates 
obtained using calibrated EMA results, this chapter first 
examines modal parameter estimation accuracy by comparing 
results obtained from an EMA on a free-free beam to those 
predicted from theory. SDM accuracy was checked by making 
various analytical modifications to the theoretical model 
for the free-free beam and comparing these to the 
corresponding theoretical results. To evaluate the 


combination of EMA and SDM, modal parameters were obtained 
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for a test structure in which the modal parameters were 
obtained only from EMA. Then using SDM, predictions of the 
modified dynamics for mass, stiffness and damping changes 
were compared to EMA results obtained from the physically 
modified structure for each modification. 

The last chapter briefly summarizes many of the 
results, assumptions and improvements as well as future 


considerations. 


2. THEORETICAL ASPECTS 
A continuous, elastic structure can be discretized into 
a system of lumped parameter elements. The motion of this 
system can be approximated by a set of simultaneous linear 
second order differential equations. In the time domain, 


these equations would be of the form: 


in ty) = ie ey eet = ee (2.1) 


where [M], [C] and [K] are the mass, damping and stiffness 
Matrices, assumed to be real and symmetric, {f} is the 
applied force vector and {x} is the resulting displacement 


vector. 
2.1 Modal Analysis Theory 


Taking the Laplace transform of equation (2.1) and 


assuming zero initial conditions gives: 


CIM] s* + [c] s + [KI] {x(s)} = {f(s)} 
which can also be written as: 


[B(s)] {x(s)} = {f(s)} (2220 


where: [B(s)] = [M] a Kelp ss 1h 


The matrix [B(s)] is called the system matrix and its 
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inverse is referred to as the transfer matrix, [H(s)]: 


[u(s)] = (e(s)]"! = oo Ris tt i) 


Equation (2.2) and equation (2.3) thus lead to: 
PHS) Eth (S) = ix(S) 0 oe C2e a) 


The Laplace variable, s, is complex and can thus be 
represented by its real and imaginary parts as Ss = 0 + jw. 
Therefore, an s-plane can be defined where o is referred to 
as the damping axis and jw as the frequency axis. The 
transter function matrix (H(s)] 1s complex’ since it is a 
function of s, and therefore its magnitude can be plotted on 
the s-plane. At particular points on the s-plane called 
poles, the magnitude of the transfer function goes to 
infinity. The poles occur as complex conjugates where each 
pair represents a particular resonance or mode of vibration. 
As a result, for an n degree-of-freedom system, there are 
n-pole pairs describing the damping and frequency of en 
modes. Each pole pair can be expressed at a particular mode 


(k) as: 


where: '*' denotes the complex conjugate, o, 1S the damping 


coefficient (negative for stable systems) and a, is the 
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natural frequency. 
it is clear that the elements of 


From equation (2.3) 
[H(s)] are rational functions in s. Therefore they may be 


written in partial fraction form as: 


ae | (2.5) 


2n 
H(s piles =f. 
k=1 9 7 PK ks] 


lan pee 


where: [a,] is a complex residue matrix containing spatial 
behavior information at each mode. 


Premultipling equation (2.5) by (s-p,)[B(s)] gives: 


a B H(s = - B 
(s - p, )(B(s)][H(s)] = (s - pp JE eee s- P, 


At a mode, s = p,. Substituting this expression into the 


previous one results in: 
[B(p,)]fa,] = 0 


It can also be shown in a similiar manner by postmultiplying 


by [B(s)] instead, that: 
Ca, J(B(p,)] = 0 


Therefore, the rows and columns of [a,] must be linear 
combinations of homogenous solution vectors {u,} such that: 


[B(p,)] {uj} = {03 
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The {u,} vectors are defined as the mode shape at the pole, 
P,. The rows and columns of [a,], must then be a scalar 
multiple of the vectors {u,}. Thus the matrix [a,] can be 


expressed as: 
+ 
he yee he (2.6) 


where 'T' denotes the transpose and A, is a complex scalar. 
Substituting equation (2.6) into equation (2.5) yields: 


af T 

n A, {u, }{u, } A* {ux}{u*} 

[H(s)] =x a oe 
k=1 Sls = ota 


(ee) 
Equation (2.6) and (2.7) show that every row and column of 
the residue matrix and thus the transfer function matrix, is 
comprised of the mode shape vector multiplied by a different 
element of itself. Therefore, only one row or column is 
needed at each mode. Equation (2.7) also shows that the 
tuansrerecunc. LON Matlixecan ubesburltCentirelyytrom the 
modal parameters (natural frequency, damping factor and mode 
shapes) and that it completely describes the modal model. 
Further reduction of the problem can be accomplished by 
computing the Laplace transform along the frequency axis 
where s = jw. This subclass of the Laplace transform called 
the Fourier transform results in a particular type of 
transfer function referred to as the frequency response 
function.’ Since the transfer function is represented 
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frequency response function if required. 
In order to measure the frequency response function, it 
is possible to rewrite equation (2.4) in matrix form. 


Assuming a two degree-of-freedom system this results in: 


pat oe GS) eames (2.8) 


ho, (s) hoo(s) )fo(s)( Xo(s)( 


where Ss = ju. 
Since only one row or column of the transfer function 
matrix is required, consideration of the first equation in 


equation (2.8) serves as an example: 
x, (s) = hy, (s)f;(s) + h)o(s)f,(s) (2.9) 


Torrorm the’ tirste row Of “the transter “function matrix, 
the response is measured only at the first point while the 
excitation is applied and measured at one point at a time to 
all points of interest. The first element in the first row 
can be found by exciting the structure at the first point 
while simultaneously measuring the response at the first 
point. Equation (2.9) then gives: 


ase i 
hy 7 () = F(a) hy, () = hy, (s) 
Ss 


= jw 
The second element in the first row can be found by exciting 
the structure at the second point while simultaneously 


measuring the response at the first point. Thus equation 
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(2.9) gives: 


x (w) 
byes iz eae) 
s = jw 
Any other row or column can be built in a similar manner and 
each element is formed from taking the Fourier transform of 
the response and dividing by the Fourier transform of the 
excitation force. 

The modal parameters can then be estimated from the 
measured frequency response function by curve fitting each 
transfer function in a row or column of the transfer 
function matrix using equation (2.5) for s = jw. This curve 
fitting process is discussed in the following chapter. 

As will be shown in the next section of this chapter, 
structural modifications are accomplished by a coordinate 
transformation to modal coordinates. Thus the following 
describes the formulation of the modal matrices derived from 
the modal parameters. 

Since only one row or column of the residue matrix [a,] 
is required, consider for example the jth column. From 


equation (2.6): 
{a, 5} = Ay {uy tuy = Ayu, 3 {uy } (25410) 
where u,,; is the jth element of {ux}. 


The ‘scalar, A,, can be expressed from equation (2.10) 
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where | | denotes the vector magnitude. 


The value of A, clearly depends on the normalization of 
the modal vector {u,}. Assuming measurements are available 
from the jth row or jth column of the transfer function 


Matrix, the modal vector may be shown as: 


{u,} = — ta, 5} (ean) 
kjj 
where ae: is the jth row, jth column element of [a,]. 


Equation (2.12) can be rearranged into: 


la, «| 
Pe we ea hes 
kK jj : Fecal 

SubSt1 cucingwinto equation (2.11) °and recognizing that up, 


is equal to unity gives the simplified expression: 


AL = aK jj 
Recognizing that A, has the dimensions, sec/mass, and from 
the assumption that the mass, stiffness and damping matrices 
are real, Richardson and Potter [26] defined the modal mass, 


damping and stiffness as: 
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For the normal mode case (zero damping), the system matrix 


is real symmetric and thus A, and p, are purely imaginary 


and the vector {u,} is real valued. Therefore: 


A, = -At 
Py = “PF 
{uy } a {uF} 


As a result, the modal mass, 


the form: 


stiffness and damping assume 
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k= ZA, (2.15) 


and the transfer matrix becomes: 


n fu, Hu}! 
(Hs) = 2 [ 4] 
Kale my saat ke 


2.2 Structural Dynamics Modification Theory 

Structural Dynamics Modification is a technique 
developed by Structural Measurement Systems Inc.. Following 
is the adaptation of this technique from Formenti and 
Welaratna [6]. 

The removal of the forcing function from equation (2.1) 


represents the free motion of a damped system expressed as: 


[M]{x} + [C]{x} + [K]{x} = {0} 


Modifications to the system can be implemented locally using 
a technique called Diakoptics or eigenvalue modification. 
This technique was initially developed by Kron [14], Simpson 
and Tabarrok [28] and Weissenburger [31] and extended by 
Pomazal [18], Hallquist [8] and Hallquist and Snyder [9]. 
Local modifications with respect to the physical system 


may be shown as: 
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where a, B and y are the values of the respective changes in 
physical units. The vector {1} locates the change in 
physical coordinates. 

This procedure allows modifications in the form of 
additions or removals of point masses, scalar dampers and/or 
Springs applied from a point to ground and scalar dampers 
and/or springs applied between two points. This same 
procedure also enables the joining together of substructures 
by scalar dampers and/or springs. 

To illustrate the local modification procedure, a 
Simple 3 DOF discrete system as shown in figure 261 can be 
used. Since the system coordinates are in a static coupling 
form, modifications can be made in physical coordinates in 


this case. 


Letting k, = 0, the following stiffness matrix is 
generated: 
Ae ey ie 
0 -k3 (k, + ky) 


In order to regenerate the original stiffness matrix, the 
local modification procedure can be applied such that k, is 
attached from point 1 =to ground. Thus, the modification 


Matrix would be: 


[ak] = k, (1341)! 
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Figure 2.1 A Three DOF System 
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where the {1} vector assumes the form: 


and: 
ky 8) 18 
k, 0} =] 0 0 0 
Omer .0 


OQ matrix regenerates the 


Adding this matrix to the k, 
original matrix. 
Supposing instead that kz = 0, then the following 


stiffness matrix would result: 


To regenerate the original stiffness matrix, the 


modification would take the form: 


[ak] = k,t1}003" 


where this time a spring is added between point 1 and point 


2 so that the {1} vector would take the form: 
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and: 
ko -k, 0 
Les 
ko {1} {1} = -k, ko 0 
0 0 0 


By adding this matrix to the kz = 0 matrix yields the 
Original matrix. | 

Since the mass, damping and stiffness changes are being 
implemented at the same point, the modified system equation 


would take the form: 
(M]{x} + [C]{x} + [K]{x} + [AM]{X} + [ac]{x} + [AK]{x} = {0} (2.16) 


Instead of using the Laplace transform as before to obtain 
the solution vector {x}, it is also possible to simply 


assume a solution of the form: 
ide 
Substitution into equation (2.16) yields: 


[a2[M] + alc] +[k] + (an? + aa + y)}1}'103 = (0) (2.17) 


Because of practical meaSurement considerations (to be 
discussed in the next chapter), a real mode model with 
proportional damping (as opposed to a complex model) is 


assumed. The mode shapes {u,} for the system are therefore 
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real and can thus be used directly in making a coordinate 
transformation to uncouple the equations of motion. As will 
be shown below, by also normalizing the mode shapes to unit 
modal mass, the uncoupled equations can be written solely in 
terms of the modal parameters. 

In order to normalize the mode shapes to unit modal 
mass, the mode shape vector {u,} is scaled such that it 


would take the form: 
{uy }/ vin 


where m, is the modal mass obtained from equation (2.13). 
The mode shape vector is now in the appropriate form to 
uncouple the equations of motion in terms of the modal 
Parameters. This is done by a coordinate transformation to 


modal coordinates using the modal matrix: 


{x} = (0]{Z} (2.18) 


where: 


{uy} {up} tu, } 
Sy ee eed aes 
ym vm vm, 


and {z} are the modal coordinates. Substitution of equation 
(2.18) into equation (2.17) and premultiplying by [e]! 
gives: | 
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Since the system was assumed to be proportionally damped and 
utilizing mode shape orthogonality, the equations of motion 


in modal coordinates become: 


he (eal area eee Seed (2.19) 


+ (an® + ea + y){o}{o} ]{z} = {0} 


where {¢} = [6]'{1} ana: 


go - Original structure damping 
coefficient 

@® - Original structure natural frequency 
Q - modified structure complex 
eigenvalue 


{z} - modified structure eigenvector 


Since the equations are uncoupled, it is possible to write 


the ith equation as: 
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From equation (2.21) it is clear that an expression for the 
kth element of the vector {z} can be determined in terms of 
any ith element. This fact leads to: 
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Substituting equation (2.22) into (2.20) yields: 
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which reduces to: 
= : mK (2.23) 
(ote See zr — ° 
(an? + BQ + y) k=l (a? + 20, 2 + ee + w°)) 


Equation (2.23) represents the modified systems 
characteristic equation. Once the eigenvalues have been 


solved from this equation, the eigenvectors can be computed 


af¢ %02 oizes3qx6 18 dant Goal> ei Fi (TS.S5) sotsatpe mere 
* rr « , es a) 2 =e 

to amass ni Banteveyes sdi’nss ‘s)} actoev sdz to tnemele aga. 

*” le “a a 


‘a9 2baol test 2id? .2nsmets Agi 


(2,3) 


=| 5 y 
(f gy’ + 7? + 4 he) ee 1) 
sof faery tos ©’) oe cys ($s e) ws + eee - ; a a 
7 ‘.— -~s* LA * Ge S05 a) AGis Geos snisuels atl 


te £ g. . 2. 
Pras pot ©) + iS + ap ‘ 


. 
‘ 


23 


from equation (2.22) to within an arbitrary constant, Aj. 
The kth element of the ith mode is then: 


oy 


262 fh eee 
k ir. 2 v) (2.24) 
(2.° + 20,25 + (a, + wy )) 

The constant A; is found by normalizing the modified mode 


shapes to unit modal mass so that: 
T 7 
{z;}'{2;} =1 


Since the modified mode shapes are in terms of modal 
coordinates, the transformation back to physical coordinates 
can be accomplished using equation (2.18). As a result, for 


the ith mode: 


ic} = [oJ{z;} 


Thus, the eigenvalues and eigenvectors of the modified 
System are determined using the original modal parameters 
and applying the assumption of a real mode model with 
proportional damping. The implementation of this procedure 


is shown in the next chapter. 
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3. PRACTICAL ASPECTS 

In the previous chapter, the frequency response 
functions, which are the basis of determining modal 
information, were obtained by dividing the Fourier 
transforms of the simultaneously measured excitation and 
response time domain signals. While theoretically this is 
Straight forward, in a practical sense there are several 
difficulties and inaccuracies which must be overcome. This 
chapter therefore begins by outlining the practical problems 
of modal analysis associated with signal analysis and 
hardware for data acquiSition, data processing and 
calibration. The final section deals with computing 
considerations for the practical aspects of SDM necessary in 
obtaining modified eigenvalues and eigenvectors from the 


modal parameters. 


3.1 Experimental Modal. Analysis 

EMA relies on broadband excitation signals in which 
most of the modes of interest are excited simultaneously. 
The various classes of these signals include transient, 
random and swept sine. In this study, the transient method 
using a hand held hammer was chosen because of its easy 
implementation in various situations. The other two methods 
rely on the use of shakers and thus are more cumbersome. The 
transient method, unfortunately, has two major 
disadvantages. If a structure is large and heavily damped, 


it becomes difficult to adequately excite it in this way. 
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Also the method is inadequate when nonlinearities are 
present since the spectrum level cannot be precisely 
controlled. The test structures in this study, however, were 
small and lightly damped with little nonlinearity so that 


these factors did not present a problem. 


3.1.1 Data Acquisition 

The experimental set-up for data acquisition using the 
transient method is shown in figure 3.1. Each component will 
be discussed below in relation to its significance in the 
Signal analysis process and its use. 

The hammer represents the basic component in the 
transient method. Hammer weight and tip hardness determine 
the excitation force frequency content. Since the structures 
considered were relatively small, a light weight hammer of 
about 200 gm was chosen. The load cell mounted in the head 
of the hammer to measure the force, had a sensitivity of 
about 2.2 mV/N. Structural frequencies below 1000 Hz were to 
be examined and therefore the tip was chosen so as to 
concentrate most of the energy in this range. As a result, a 
hard rubber tip was used, which resulted in an excitation 
force spectrum that was flat (+3 dB) to about 1300 Hz. While 
the input force was being measured by the hammer~-mounted 
load cell, the response was simultaneously measured by an 
accelerometer. An initial accelerometer used came with the 
hammer kit, but was found to cause accelerometer loading 


errors. The accelerometer's 26 gm weight was substantial 
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enough to affect the dynamics of the test structures. Thus a 
much lighter accelerometer (2 gm) was used. While this lower 
mass gave a small accelerometer loading error, it also 
reduced the output sensitivity from 10 mV/g to 2.37 mV/g. 
Since analysis or processing of the excitation and 
response time waveforms produced by the hammer and 
accelerometer requires digital computation, the signals must 
be digitized and limited to a finite period. The discrete 
waveforms thus have a finite resolution (At) based on the 
number of samples (N) and the period (T). This time domain 


relationship can be expressed as: 


T = Nat 


In order to transform the digitized waveforms to the 
frequency domain, the Discrete Fourier Transform (DFT) is 
used. The algorithm which accomplishes this task, is called 
the Fast Fourier Transform (FFT). The resulting frequency 
waveforms are complex quantities and are thus made up of two 
parts, namely the magnitude and phase. Therefore, N time 
domain samples will yield N/2 complex frequency domain 


samples. This frequency domain relationship may be shown as: 


where F,, is the maximum value of the frequency in the 


frequency range of interest and Af is the frequency 
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resolution according to the relation: 


where F, is the sampling frequency. N is a fixed value in 
the FFT algorithm (typically N = 1024). Therefore to extend 
the frequency range of interest, the sampling frequency is 
increased and thus the time period is shortened. Note then 
that this also results in lower frequency resolution. 

From the previous equation, it can be seen that the 
sampling frequency, F,, must be twice the maximum frequency. 
This relationship is referred to as Shannon's Sampling 
Theorem. If F, is less than twice, a problem called aliasing 
can occur. This problem refers to the fact that signals of 
frequency greater than the difference between F, and F, will 
overlap or alias in the frequency range of interest. 
However, by setting F, at least twice F,, the difference 
would then cause the alias signals to fall above the F,, 
limit. Since aliasing can also occur due to stray signals 
from surrounding sources, a low pass filter is used to 
ensure that the frequency range is limited. Because of the 
Finvce tiiter roll-oLrt (a Characterrstic of all-friters) =the 


cut-off set at F, does not ensure that all higher 
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frequencies are completely rejected. Signals occurring in 
the roll-off band might not be entirely attenuated by the 
time they alias in the frequency range of interest. To 
acCcountmlorethias See, the sampling frequency is uSually 
raised to 2.56 times the maximum frequency. The filter used 
ineehasestudyuhad asroid-off rateso£ 96udB/octavessThnus, if 
F, = 2.56F,, then alias signals would be attenuated by 
approximately 62 dB. This is more than enough since the 
entire range considered in this study was 60 dB. 

Digitizing is done by an analog to digital converter 
(ADC). The type used was a twelve-bit ADC with an effective 
Gangenof ’72°dBvor +57 volts. Astavnresult; /tocmaximize the 
Signal-to-noise ratio, an amplifier was needed to boost the 
Signals thereby filling up the required range. It was also 
necessary on the other hand to ensure that the signal did 
notesaturatesan?the range or) -elipping'’ .ofethessignali and 
thus distortion occurs. Another consideration was that for 
points on the structure where there was a lot of response, 
it frequently occurred that after the initial hit, the 
structure would spring back into the hammer before it could 
be moved away resulting in a multiple impact. The impacts 
after the initial hit are typically of such low level at 
particular frequencies that the signal-to-noise ratio is 
very much reduced at those points in the frequency spectrum. 
In order to ensure that the signals properly filled the 
required range and that multiple impacts had not occurred, 


the digitized time domain signals were plotted by computer 
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graphics for viewing. Acceptable signals were then stored on 


floppy disc for future processing. 


3.1.2 Data Processing 

Since the time domain data is a limited sample taken 
from a continuous waveform, the data thus represents the 
viewing of the continuum through a rectangular window. This 
corresponds to multiplying the time domain data by a 
rectangular window function. Multiplication of time domain 
data by a window function is equivalent to convolution of 
the Fourier transforms of window and data in the frequency 
domain. Distortion of the transformed waveforms will thus 
result unless the time signals are periodic (waveforms form 
integer multiples of cycles in the finite time period) or 
decay completely in the time period considered. Added 
distortion can also occur since random noise in the 
measurements is not periodic. This distortion is referred to 
as leakage since signal power 'leaks' out over the frequency 
spectrum when transformed. Various types of window functions 
can be used to minimize the leakage depending on the type of 
waveform examined. Windows also have the effect of reducing 
the amount of noise ina signal. 

The time domain window transformed to the frequency 
domain assumes the form of a large main lobe followed by a 
series of smaller side lobes that occur over the rest of the 
frequency range. Signal distortion is minimized when the 


main lobe width and side lobes magnitude are minimized. An 
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unfortunate conflicting property of windows is that by 
decreasing the main lobe width, the side lobes magnitude 
increases and there are less noise reduction effects. These 
properties are proportional to an increasing time domain 
window width. Thus compromises must be made in choosing the 
proper window. However, a standard requirement of windows is 
that to minimize leakage, the window should go to zero at 
the ends of the time record thereby forcing the signal to be 
periodic. For example, a common window used on sine waves or 
random signals is called the Hanning window. This window 
resembles a Bell curve which tapers to zero at both ends. It 


can be expressed as a function of the form: 


w(t) = 0.5 [1 - cos(2nt/T)] Op<atac) 

When the impact method is used, windows take on an even 
greater importance. An important technique used as well in 
reducing the noise is to take the power spectrum and 
cross-spectrum of signals at the same points and average 
them. These are then used to form the frequency response 
function (to be discussed below). Since the whole point of 
the impact method is to save time, less averages are taken 
and thus noise reduction by windowing is very important. 

Although the forcing signal for the impact method is 
self-windowing since it dies out in the time record, stray 
Signals can occur due to hammer movement after the impact. 


One window which yields reasonable results in this case, is 
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the squared cosine which drops from unity to zero in 1/16 of 


the time record. This window can be expressed as: 
w(t) = cos@(8nt/T) 0<t<T/I6 


For the case of the response signal in an impact 
excitation, an exponential window is often used. If the 
Signal decays rapidly, then much of the sampling time period 
will contain noise. On the other hand, if the structure is 
lightly damped, then there may be very little decay of the 
Signal and thus there may be no periodicity. This effect 
will then produce leakage. Both these problems are typically 
solved by multiplying the response signal by an exponential 
function which decays exponentially from unity to 0.05 in 
the time period. Since the important information occurs at 
the beginning of the time period, it is weighted more 


heavily. The exponential window is represented as: 
w(t) = ert ORs teal 


where: 
hee gn(0.05) 
i 
This particular window has the effect of adding additional 
damping to each mode. However, this damping is a known 
amount, b. Thus, once the modes have been identified, the 


damping can be simply subtracted from them to obtain the 


actual damping. One unfortunate consequence of this window 


ssgml nk nt feqcle Seeeges2 ak7 Yo saeco ofS 

edd 2f .Beey nedio' 2: gebaiw [siiastogte ae ,Hol 
hoitag sais prilqene-sd) 3q fovm asts »ylbige? zywosd f | 
ef ewitsursea ad> ti fares aacJ9q aria co .Seian aissao3 p 
gd3 Fo yaoSb S1sS1L Visv ea tyem svois vets \begash Yla 
Sestio ‘e227 SutibiSai iy on a: Yam axed? aode Sa6.% sa ies 


; 7 _, 
gilssiqvue 215 awalcarg*seetls aycl -spxtasl suubedg ned? Ifiw - 
feizmenegxs as yo Deng) 2 eanogeas saz ontylgesiog ya havioe 
7 


> 
ol 20.0 62 vial mevi @flsidnanoges overs’ doinw mezgon ‘ 


- 


$6 @3usse hertessctnl Jnesicdmi eat eante .Bborieg anita? 5 

droft batdpisw es ts “Soiveq aris edy do. enintiped add 
> ¢ » ~ . : cr a “ 

rag Saroeseicet Jai -wohniw Lalinsdigks et? -yLiveed 


Post 30 419 © fae 


fametzipis pniibe to 259823 ede ee WebAES 


_ 


o 


33 


occurs if the system is already heavily damped. This 
additional damping may cause the modes to become even more 
closely coupled and thus harder to identify. In these cases 
Or in general where damping is to be determined more 
accurately, a Zoom transform can be used. The Zoom transform 
increases the resolution by concentrating on a narrow 
frequency band instead of the entire baseband. This 
technique is discussed in Ramsey [21] and was not used in 
the present study. 

It is important to note that with both the forcing and 
response windows, there was no need to force a zero 
amplitude at the beginning of the time period. Both are 
inherently zero as the structure has not yet been excited. 

The determination of the frequency response function 


(FRF) can be determined from its basic definition: 


However, to reduce the inherent noise effects, it is 
actually done by using averages of power and 
cross-spectrums. To understand how this is done, consider 


the response spectrum as: 
ee ee eS 


where Ss, is the noise Fourier spectrum. While the effect of 


the noise can be substantially reduced by averaging the 
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spectrums, since they are complex, they do not readily lend 
themselves to this averaging process. It is easier to use 
the power spectrums which for the response and forcing are 


defined respectively as: 


G.. = S.S# Gee 5.5% 


Taking the square root of the ratio of the above quantities 
would yield the FRF. However, in this form there is no phase 
information and the noise will not average out. The 
cross-power spectrum can be used for this purpose. Tt can be 


defined as: 


Averaging in the presence of noise yields: 
Grp 7 SSE = (SH + SSE = Geet + Gog 


and thus: 
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where the second term numerator is the noise cross~-power 
Spectrum. Since the noise and forcing signals are 
uncorrelated, the noise cross-power spectrum will average to 
zero and thus a closer approximation to the true FRF can be 
determined. 

In order to quantify the amount of noise in a FRF 
measurement, the coherence function is computed. This 
coherence function can be defined as the ratio of the 
response power caused by the applied force to the measured 
response power.The measured response power can be expressed 
as: 
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The response power caused by the applied force is: 
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where averaging will give: 
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Therefore, the coherence can be expressed as: 


As a result, if the coherence is one, then the response is 
due solely to the forcing and thus no noise is present. All 
values of less than one suggest noise contamination 


quantified by the particular value. 


3.2 Modal Parameter Estimation 

Once the FRF's have been formed, accurate estimation of 
the modal parameters is done by curve fitting. Each FRF 
formed at each point on the structure represents the sum of 
the contribution of each mode in the frequency range of 
interest. If there is a large amount of frequency separation 
and little modal damping, then the effect of one mode on 
another is minimal and a curve fit that fits one mode at a 
time is adequate. The assumption then is that each mode can 
be represented by a single degree-of-freedom model. However, 
often the situation is such that the modes have a marked 
effect on one another and there is thus strong modal 
coupling.. In this situation, the modal parameters must be 


identified simultaneously. Single mode and multiple mode 
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methods are discussed in some detail in Richardson [24]. A 
very good multiple mode curve fitting algorithm is described 
in some detail in Richardson and Formenti [25]. The method 
used in this study for modal parameter estimation was an 
algorithm developed by Fyfe [7] based on a complex curve 
fitting scheme devised by Levy [15]. This method curve fits 
the experimental data to the analytical function as 
POEmMUlatedsineequacvonm¢c.5)) [Ores = ja.) Ine wit.) 16nd 
multiple mode one in which minimization of error uses a 
least squares technique. This method is thoroughly discussed 
in the reference by Fyfe. 

It should be noted that usually modal analysis revolves 
around the use of specialized equipment such as the Fourier 
analyzer to collect and process the data. For advanced 
analysis like multiple mode rect e cana a computer is 
also required. This study peeeeed the Fourier analyzer and 
used a complete mini-computer software package developed by 


Fyfe to collect, process and curve fit the data. 


3.3 Calibration 

The traditional technique of calibration involved the 
determination of the individual sensitivites of the 
accelerometer and load cell by comparison with a standard. 
This method, however, has several problems associated with 
it. These include inconvenient on-site calibration, 
accumulation of errors due to individual calibration and the 


fact that amplitude is usually measured at only a few 
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frequencies with phase being entirely ignored. In addition, 
there is typically a lack of regard for auxillary equipment 
such as amplifiers and filters. 

For the impact technique, the traditional calibration 
can lead to entirely erroneous results. The force picked up 
by the load cell is not equal to the force exerted on the 
structure due to an effective mass between the load cell and 
the structure. This effective mass being due to the impact 
tip mass and its material properties. Another technique 
which circumvents many of these problems is called ratio 
calibration. The combination of hammer, load cell, impact 
tip and accelerometer used to measure the dynamics of the 
structure considered can be calibrated simultaneously. This 
task is accomplished uSing a very simple experimental 
Set-up. 

A relatively large block of known mass is hung from a 
long suspension. This set-up can be made to approach a 
Single degree-of-freedom system if the suspension is long 
enough and if there is little motion through the block's 
center of gravity. Thus by mounting an accelerometer on one 
end of the block and impacting it through its center of 
Gravity at. thesopposite. end, the) force ican) be,expressed 


through Newton's law as: 


F= mA 
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This relation thus represents the inertance or expected FRF 
of the calibration block and is a constant (inverse of the 
block mass). A calibration function can be defined then as 


the ratio of the expected FRF to the measured FRF expressed 


as: 
H = ve = ] 
c A/F, mH) 
wheres: 
Ho - calibration function 
Hy - measured FRF of calibration block 
Mm. —~smassS of calibration olock 


Assuming linearity for an actual structure being tested, its 
measured FRF can be calibrated by multiplying it by the 
Calibration function to obtain the actual FRE. This 


relationship can be simply shown as: 


The actual FRF will be in physical inertance units depending 
on the physical units used for the calibration block. To 
convert to compliance units, it is only necessary to divide 


through by, (ja) *. 


on aa abzovat) suetameo 88 
an aada peaibek ‘od! 6S noisy 


bashasexs at5 bawased ott oF ae een srit Yo ¢ 


awed ial te : 


noiganw2 noreeedtiss - # 


%56ld- “el ¥stdliss. Ia FAY Be=uwacga - 7H : 
f o i a a - 
A20id RO. TSIGLisS 1S Eee. 


asi ,betees oniad SxuigouTse Leuds6 fs- tol yo reese 
adj vo at eoivloisium yd betewdiies sd ca6 SRS Ber 
—— © 


2iAT . 4%. igutas ois nigdadn os gosvonus ai 3e3¢ tiss 


ten Avene ~igmic-at thS ciieneisels 
Hig gt 


nnbasdsa esi sanes.er. pas ares, ac. is iw 
‘ 


Oo? .toold “eiserdtle> S43 4of baeu aetay — 


 sbivib-oor yssagezen Yee @igh ,aginw sonal § 


Sth 


40 


3.4 Structural Dynamics Modification 

Structures whose mode(s) possess proportional or very 
light damping can usually be described in terms of a real or 
normal mode model. At a particular mode, this model is 
characterized by a mode shape which attains maximum or 
minimum deflections simultaneously. As a result, different 
points on the structure are either in phase or 180° out of 
phase and thus nodal points are stationary. In actuality, 
all structures possess damping where the modes are complex 
to some degree. This implies that a mode shape will not 
reach maximums or minimums simultaneously and thus nodal 
points 'travel' i.e. the phases will assume values other 
SnaneG mor 160. 

However, Since a majority of the mechanical structures 
dealt with possess relatively little damping, it is not very 
often when truly complex modes are encountered. Even when 
they are, they are almost impossible to adequately identify 
as such due to practical measurement problems. These include 
the inherent problems of noise, transducer and 
instrumentation phase match characteristics, leakage and 
errors in parameter estimation. Formenti and Welaratna [6] 
examined several of these problems and their effects using 
Simulation. They have found through experience that the 
errors in assuming a real mode model are much less than 
those occurring when trying to estimate complex modes. As a 
result, the SDM technique utilizes a real mode model. 


Damping is taken into account through the assumption of 
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proportionality as shown in section 2.2. 

The curve yiitter in this study (section 3.2) results in 
complex frequencies and mode shapes. Therefore, the mode 
shapes can be approximated as real by considering only the 
magnitude of the residues and determining their sign by 
rounding off the phase. Because of the assumption of 
proportional damping, the complex frequencies can be 
sohadnel as complex. When the modal parameters have been 
converted, a real mode model with proportional damping can 
BbeUconstructed. 

Once the physical modification value and location are 
decided and the mode shapes have been normalized to unit 
modal mass, then the modified systems damping and natural 
frequencies or complex eigenvalues (2) can be calculated 


from equation (2.23). Rearranged in terms of a function, it 


becomes: 
4. 1 
AYA ar 2 Se anes A OPO Re a IO 
k=1 (Q + 24,2 + (a + wy )) {aQ” + BQ + y) 
where: 


{9} = [0]"¢1} 
Thus at a point on the structure where there is a mass 
modification and/or stiffness and damping changes applied to 
ground, the vector {¢} is the row of the modal matrix [€] 
COLreSpDoOndIng toOuthat particular point. or DOF! For danping 
and stiffness applied between two points or DOF Js, the 


vector {¢} is the difference between the two correspondina 
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rows in the modal matrix. 

Equation (2.23) represents the modified systems 
characteristic equation. It is possible to expand equation 
(2ia2a:) tintoeascharacteristicapolynomial. Thedsolutionscould 
then be accomplished using any of the numerous polynomial 
algorithms available. However, for systems greater than 1 
DOF, the expansion into the polynomial becomes an 
increasingly difficult task. Fortunately, the roots can 
Still be found from the original form by using the 
Newton-Raphson method. This method is iterative where the 


iterative equation is formulated as: 


grt] - Qh 2 F(a") 
F'(Q") 


y \J 


where 'r' denotes’ the iteration and F’ is the functional 


derivative of equation (2.23) expressed as: 


. 2 
n =2d mn et, ) 
, k k (202 + 8B) 
F(a) = Sh = [ r —_*——_*, 5] - ec 


J 
k=1(9° + 20,2 + (o,° + 02) (an? + ga + y)é 


Once the complex frequencies have been determined, it 
is a straight forward process to obtain the eigenvectors or 
mode shapes by direct substitution into equation (2.24). 
Transformation of the modified mode shapes back to physical 
coordinates is then accomplished by uSing equation (2.13). 

The Newton-Raphson method was both easy to program and 
converged to complex frequencies rapidly. However, for very 
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modification values, it was found that convergence to a 
modified structure frequency did not occur. As a result, it 
was deemed a valid assumption that the modal frequency was 
not affected and thus the previous unmodified value was 
retained. 

Initially the Newton-Raphson method was coupled with 
the Incremental-Search method so as to make the root finding 
process automatic. It was found to be much faster however, 
to simply use an interactive approach whereby initial values 


were manually input. 
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4. RESULTS AND DISCUSSION 
In the previous chapters, the theoretical and practical 
aspects of EMA and SDM have been individually examined. This 
chapter follows the next logical steps of first separately 
verifying these techniques through an example of a well 
known structure (a free-free beam) and then evaluating the 
combination of the two techniques by application to a more 


complex structure. 


4.1 Calibrated EMA Validation 

As a check on the modal analysis technique, a free-free 
beam was experimentally tested and the calibrated results 
obtained were compared to theoretical results. The 
calibration procedure followed that as outlined in section 
3.3. The calibration mass used was an arbitrary block of 
steelewith dimensions! 2eanu x 2) ine x 5-in (50. 8mm) x) 502/65 mm 
x 127 mm) resulting in a mass of 0.0146 lb-sec*/in (2.55 
kg). Initial theoretical calculations for the free-free beam 
placed the first three modes at just under 1000 Hz and as a 
result, the calibration FRF was determined with a cut-off at 
WOOVHZe YI tei stimportant™ to: notemthatesinces two FRE sowere 
multiplied together (section 3.3), they must have compatible 
frequency resolutions. 

The experimental set-up for the free-free beam 
consisted of a steel beam with dimensions 3/8 in x 3/2 in x 
22 in (9.53 mm x 38.1 mm x 559 mm) suspended by elastic 


bands. Eleven equal locations from one end to the other were 
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marked Off and tested for lateral vibration up to. 1000 Hz. 
The same accelerometer as that used in calibration was 
mounted at the first point at the end of the beam in the 
vertical direction and the same hammer was used to impact 
the beam at all eleven locations. 

The first test point coherence and calibrated inertance 
FRF is shown in figure 4.1. The entire plot actually 
consists of the three separate plots of coherence and the 
FRF magnitude and phase. The dashed line corresponds to the 
curve fit results. It can be seen from the figure that there 
are three well defined modes with very little damping. Since 
the coherence is essentially unity at all three peaks and 
the curve fit conforms to them very well, this measurement 
Should be expected to accurately portray the dynamics at 
this point. Table 4.1 shows the excellent agreement between 
the experimental (averaged over all the test points) and 
theoretical frequencies. 
Since SDM uses mode shapes normalized to unit modal mass 
(section 2.2), and since this directly corresponds to the 
degree of modal participation as well as calibration and EMA 
accuracy, both the calibrated test and theoretical mode 
shapes were compared in this form as shown in table 4.2. The 
table shows that the first and second mode shapes are very 
close to their theoretical counterparts while the third mode 
shape is not as good. It was found that in most of the 
measurements, the coherence was not as good in the vicinity 


of the third mode. This fact is especially evident at test 


aia ita masd Virago silt 3 3 Bia ait aig 4a 6 is a 
toaget os Sseu nO remeaBs 8 nee Sig ne noitsers& Lod 7s z 


-2na7assol revels Lis 36 arr 


soneaioeni beaescilss ona oaertites. Shing 2282. sexid edt 
yilguges joig ovtiae att ott eqnpit ai awods et) 


a 
= 


eds bis spnstenoo to ade te sigiaqse seta odd Io etme 


$13 ¢2 ebnsaesz209 ens, Hereat Sa7 Sead Dos stuainpss: 
siedt =e “tleii siz maxi aeSe “sd oad 3¢ -etivess $33 
egnie .pniameb efagtt yrewids tw sagem bettie ' law easdz . 

Snes Geeeq caudt Meta gdinu yitetsnsees et 2one26ae3. 


inemesveeem 22h7).lisw~ yvisv mais ci zwieinos 36? Sree 


& enimenuh sit ysicitG Yfasevests, 03  Secsequs acd 


» . Bei oreimend beviseae 

seem Lahom sia of basiiamion esqaie .shorm. spac Mee 

si3 .o7 Hadedd-ficd ¥ieHbst® wbx Mot brs (8.0 ee 
aus Sos netdssdiiss 25) [lew as naitsoizisieq i shom- ae 
aton Laoideyosid bas tee Bedsidilss sda dso ye 

sdT .sit-sidad- nl nwotle a6 m103 ecdt nt bstagnoe sisew 
Use” 3x5. sstrade abom, baovse hag 22273 ant Sens 
show Sains ads -elinw ees an seria 
ad3 Yo seem af gsdayba hdl 
qsietoty mee a “por ap 


‘ 7s 
= nae AG te: Ries ae Pes 19 nt 3% 
> a a -~ aoe > a b= ° 


ey ne: 
$378" ale 7 
7 


ean 
See bd 


-DAMPING (%) FREQUENCY (Hz) 
MODE THEORY TEST 9 THEORY TES 


1 0.0 0.06 19820 159.4 
2 0.0 0.06 436.8 438.8 
8) 0.0 0.04 plese! 858.8 


Table 4.1 Comparison Between Theoretical and Experimental 
Frequencies for the Free-Free Beam 


LEST MODE 1 MODE 2 MODE 3 
POINT THEORY ee THEORY EST THEORY TEST 


] Zine 22 st ral RAY rae, Zila Lome 
2 11.4 Tera 4.8 4.6 -11.0 Seely, 
3 Zl 203 - 8.4 = 85 =|'3'.9 -15.8 
4 = 5/0 = 9.6 -14.0 -13.7 - 8.4 - 8.0 
5 -11.0 -10.8 =) 0x2 -10.3 6.9 TESS 
6 -12.9 -12.6 0.0 Or 1.1 1852 
d, <1 720 =10a7 Wee 102 6.9 oad 
8 5.0 = 0150 14.0 Tn - 8.4 -10.8 
S) Ze ee 8.4 8.1] -13.5 -16.0 
10 11.4 Tl - 4.8 - 4.7 -11.0 - 1.6 
am 2bee (false -21.2 -21.6 rans 205 


Table 4.2 Comparison Between Theoretical and Experimental] 
Mode Shapes for the Free-Free Beam 


‘ 


3 Bad! 


oe 


ES rare 
fajnemtisgx3 bie rgthionet was ncetragnad 1.4 staeT 
mso8 43%) wee 361 2atonaupen4 


3008 3 30m 130 é 
T2347 YROTNT Teast... TAOSHT T23T Wage 


5.58 5: fs f-ss ee ©. 56 a | 
‘ai Or FT. b.% “th we ait 
B,cl- ¢tf- ‘3 - b3o - ae Ts 
0. = a (eG ls OcdT- 2.2 - &.8 = 
ma oa pare SogP. @ Of- 0. Tf 
>. 8 yet bg 6.0 >. $i- e.$f- 
‘8 G3 50 conf eu) 0, ifs 
§, bt. hobo Cae Oa. a2.- g.2 = 
Gidi- Git j= ye ie nit . 
B.f = QD ii- \.&- Bek - Le P 
Bas S15 a, 1S- ome e fs S48 


| 
. 
n 
: 
m4 


tagnamnfaeexa bre Veta aendes nod Praca sh aad 
masa sav3-s59 grt vot 2aqend Shah : 


47 


FREQUENCY (Hz) 


Figure 4.1 FRF and Coherence for the Free-Free Beam at Test 
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Figure 4.2 FRF and Coherence for the Free-Free Beam at Test 
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point 2 which gave the poorest mode shape estimate. The 
coherence and FRF for test point 2 is shown in figure 4.2. 
It is evident that while the EMA technique and calibration 
can be quite accurate, noise can greatly affect this 


accuracy as evidenced by the coherence in figure 4.2. 


4.2 SDM Validation 

In order to check the accuracy of the SDM method, the 
theoretical modal parameters for the free-free beam 
described in the previous section, were used as input to the 
SDM program. The theoretical modal model used was described 
by Six modes including the two rigid body modes. This model 
is shown in tables 4.3 and 4.4 where the mode shapes have 
again been normalized to unit modal mass. This model was 
modified using SDM and compared to the corresponding 
theoretical solution. The modifications consisted of the 
separate addition of a pin, a stiffness,of 10,000 lb/in 
©1275 10ch N/m) and a mass. otm! 1b \(0k454.kq)e tol pointed: 
The pinned modification was modelled in the SDM technique by 
adding an infinite stiffness. This simply results in the 
left handisidesotmequation ((2.23)..g0ing sto. zere.6 Table 425 
represents the theoretical non-rigid body frequencies in 
comparison to the predicted using a varying number of input 
modes for the three separate modifications. In each case, 
the relative modification was large in order to adequately 
test SDM accuracy. As was expected, the accuracy diminishes 


for predictions at higher modes (since the model is 
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Table 4.3 Theoretical Free-Free Beam Frequencies 
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Table 4.4 


MODE 2 MODE 3 MODE 4 MODE 5 


-18.4 Ze epee AleZ 
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0.0 -12.9 00 Tore 
cia =11.0 1052 6.9 
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Theoretical Free-Free Beam Mode Shapes 
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truncated) or if fewer modes were used in the prediction. 
However, it 1S important to note that the accuracy did not 
diminish exponentially but much more gradually. 

. To determine mode shape prediction accuracy as well, a 
fourth change was made. The full six mode free-free beam 
model was used to predict the effect of two local 
modifications which pinned both ends of the beam. The 
comparison in terms of both frequency and mode shape 
(normalized to unit modal mass) with the theoretical is 
shown in tables 4.6 and 4.7. It can be seen that the 
prediction iS quite good overall, but accuracy again 


diminishes for higher frequencies. 


4.3 Unmodified Test Structure 

The structure chosen to be tested was built to be 
Similiar to that used by Herbert and Kientzy [11] and Ramsey 
ana Firmin [23]. “ApySchematic diagram ‘of this structure is 
given in figure 4.3. It consists of three aluminum plates 
bolted together along the adjoining edges. In order to 
determine the mode shapes accurately, a grid comprised of 
seventeen test points was established on the top plate. No 
points were used for the back plate since the motion could 
be deduced from top plate information. Only two points were 
deemed necessary for the bottom plate since it was large and 
heavy relative to the other two and thus little response was 
expected. The total number of test points then was nineteen, 
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FREQUENCY (Hz) 


MODE THEORY SDM 
1 69-9 7 Oiend 
2 27950 281.8 
3 629.1 645.7 
4 118.5 1147.7 


Table 4.6 Theoretical and Predicted Frequencies 
for the Pinned-Pinned Configuration 


TES MODE 1 MODE 2 MODE 3 MODE 4 
POINT THEORY SDM _ THEORY SDM THEORY SDM THEORY SDM 


1 0.0 0.0 0.0 0.0 0-000 0.0 
a 4.6 4.7 8.8 ele aR kiss ans) 
3 8.8 B20 11438 14-39) \4e3° 3147.3 Gao 
4 TO al2e2 8148 14.3 4.0 5.7 = 98:8 
5 14.3 14.3 8.8 Sif) =" G400=50.0) =)423 
6 1550.5 15.0 0.0 0.0 -15.0 -16.0 0.0 
7 Vase a2 8 - 8.8 =P O58 SSI 7-2 O56) 143 
8 Waly 222 9-14-38 =14.4 G26 2528 8.8 
“) 8.8 8:8 -14.3 -14.3 1453 14.3 -(6-38 
10 4.6 nO et Oc O ue -b OO MNP lems Vool modi 
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Table 4.7. Theoretical and Predicted Mode Shapes for the 
Pinned-Pinned Configuration 
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of 57. Since many of the points were internal to the 
structure, this total was reduced to 30. In order to 
Simulate a free-free condition, the structure was placed on 
a block of foam rubber for all measurements. 

As before, the modal information was obtained using the 
impact method with an accelerometer fixed in the z direction 
at point 1 (thus point 1z). Preliminary tests displayed the 
existence of five distinct vibrational modes under 400 Hz. 
Therefore, with a sampling frequency of 1024 Hz, the 
frequency resolution was a very acceptable 1 Hz. It was 
possible to impact the structure at all points in the 
perpendicular or z direction. To obtain information in the 
plane of the structure (x and y directions), selected points 
were impacted parallel to the surface. Since the 
accelerometer was fixed at point 1z, and impacted at all 
others, the first row of the transfer function matrix was 
measured. 

The FRF information which was curve fitted for point 1z 
up to 400 Hz (figure 4.4) shows the five distinct 
vibrational modes as well as the occurrence of two rigid 
body modes at the lower end of the frequency scale. The 
coherence at this point and for most of the others in the z 
direction were very good. It was however, typically much 
poorer for measurements in the x and y directions. This 
effect can most likely be attributed to low response levels 
(and thus lower signal-to-noise ratios) in those directions 


and perhaps transverse sensitivity effects of the 
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accelerometer. 

The modal parameters obtained in the z direction are 
shown in tables 4.9-4.14 where the mode shapes were again 
normalized to unit modal mass. Since damping and natural 
frequency are global properties of a structure they have 
been averaged and presented as shown. The geometrical 
representation of the mode shapes for the five modes are 
shown in figure 4.5. The first mode at 48.2 Hz represents an 
almost purely flapping mode of the top plate (first bending) 
whereas the second mode at 171.3 Hz displays almost pure 
LOrsvonmeor the top plate around the * axis (Glirst torsion). 
Thestnicdvat@zi/.9) zis scomprised or agcombinatpron. of top 
plate torsion around the x axis and back plate torsion or 
twisting around the z axis. The fourth and fifth modes at 
290.8 Hz and 395.8 Hz respectively, are essentially 
‘rippling' modes of the top plate (second bending in and out 
of phase). 

An examination of the mode shapes indicates that point 
3 undergoes relatively large excursions in the first three 
modes. As a result, this point was chosen for modifications 
as it appears that this would cause large dynamic changes 


and thereby give a thorough test of the SDM concept. 


4.4 Effects of Mass Modification 
The simplest modification to implement consisted of an 
aggre uoneOf y0.25 Lb) (mass. Of 0. 11 ka) of lead at vpoint 3. 


This configuration was tested using EMA at the same test 
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POINES aS in theoriginal structure. The MRF, curve fit and 
conerencé fare shown again for point wa in faqure 4.6. The 
most obvious effect of the mass was an overall reduction of 
the natural frequencies. 

In order to model this change using SDM, the mass was 
added forall” 3)D0OP (x,y and zidinections) at pointe 3 sor that 
in affect three local modifications were performed. The 
modal frequency comparison between the test and the 
predicted results are shown in table 4.9. The predicted 
effect on the mode shapes in the z direction is summarized 
in tables 4.10-4.14. It can be seen from table 4.9 that the 
prediction is good when considering the magnitude of the 
changes in frequency. The second and third mode frequencies 
show the greatest discrepancies most likely due to a 
combination of several factors. The first factor relates to 
the fact that these modes had the largest contributions from 
all three directions with the result that the motion is more 
complex and thus more dependent on three measurements 
instead of one. Secondly, the measurements in the x and y 
directions were more prone to noise effects as explained 
above and thus poorer estimates of modal parameters were 
used in the prediction. The last tactor,; which) 1s) probably 
the most important, is that rotational inertia of the mass 
1S not taken into account. |}Gomplete motion atva point 
involves not only the translational motion in the x,y and z 
directions but also rotations around each axis. Thus, to 


accurately predict the motion, rotational DOF's should be 
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Figure 4.6 FRF and Coherence for the Mass Modified Structure 
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added. This is especially important for beam structures, 
such as the test structure considered here, where rotational 
effects are most pronounced. Unfortunately, rotations are 
difficult to measure and thus there are no commercially 
available transducers. However, rotations can be estimated 
uSing several translational transducers. Discussion of 
rotational DOF and their estimation has been done in 
references [2] and [29]. 

A more complete way of displaying the entire prediction 
in terms of frequencies and mode shapes is to synthesize an 
element of the modified FRF matrix and compare it to the 
actual FRF measured. This synthesis is in the form of a plot 
of the magnitude for point 1z and is shown in figure 4.7 for 
the mass change. The solid line represents the magnitude of 
the measured FRF and the dashed line displays the predicted. 
As 1S seen, the magnitudes are also quite close which 
implies an accurate prediction for the mode shape as well. 
Again, the largest deviations appear to occur at the second 


and third modes due to the reasons previously discussed. 


4.5 Effects of Stiffness Modification 

The second type of change considered was the addition 
of stiffness. Once again, in order to provide large relative 
changes in the test results, the change included attachment 
at point 3. A brass rod 1/16in (1.6mm) in diameter was 
attached between points 3 and 19. Using simple beam theory 
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61 
three directions was calculated from: 


k, ; = Kok, 7 = 
L 

where Amisethe cross-sectional area, L is the length, I is 
the area moment of inertia and E is the Young's modulus. 
This rod gave calculated values of k, of 10,500 lb/in 
(1.84 (10°) N/m)yand ky and ky of 0.4 1b/in (70 N/m). The 
total weight of the stiffening rod addition was 0.016 lb 
(mass of 0.0073 kg) and this was evenly divided between 
PoInts 3%and 19..-The FRE, curve fit and coherence for the 
stiffness addition for point 1z is shown in figure 4.8. The 
main effect of the stiffness was to raise the natural 
frequencies so that only four modes occurred in the 
Ereauency range of interest. As seen in table 4.9, the 
lowest mode occurred at approximately 106 Hz whereas in the 
Original structure it was at 48 Hz. 

To model the change using SDM, the stiffnesses were 
added in the respective directions between the two points 
and also half the total mass was added to each of the points 
3 and 19. However, it was found that the small amount of 
mass and the x and y stiffnesses caused negligible changes 
in the structure's dynamics. The longitudinal) or 2 direction 
stiffness though, resulted in large changes in frequency as 
shown in table 4.9. The changes in the mode shapes can again 
be seen in tables 4.10-4.14. The tested and predicted 


results from table 4.9 agree very well. This good agreement 
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is due to the fact that essentially only one local 
modification was needed, that rotational DOF's could be 
safely ignored and that very good original modal parameters 
(z direction only) were used in the prediction. The 
Synthesis plot for the stiffness change is presented in 
figure 4.9 for point 1z. The mode shape can be seen to be 
very close as well. 

It should be noted that since five modes were utilized 
in the modal model, that a fifth mode was also predicted. It 
was found that the predicted fifth mode occurred at 1019 Hz. 
Although the actual structure was not tested for the fifth 
mode (it exists outside the frequency range of interest), it 
seems highly unlikely that the fifth mode would jump from 
395.8 Hz to over 1000 Hz. The explanation for this behavior 
is that the modal model used is a truncated one which was 
built around intormation only from D€uto 400 Hz] Theretore, 
it would not be expected to predict the dynamic effects very 
far outside the scope of this model. 

Comparable results for similiar mass and stiffness 
modifications to a similiar structure can be found in 


references [11] “and [23]. 


4.6 Effects of Damping Modification 

The final modification considered was the introduction 
of increased damping. This change was expected to yield the 
least predictable results for a variety of reasons. These 


include the fact that proportional linear viscous damping is 
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assumed and also that a real mode model is being used to 
model a complex mode phenomena. As well, it is known that 
material properties such as shear modulus and damping 
coefficient are quite dependent on frequency, temperature, 
dynamic load and static preload. 

The physical modification consisted of a damper 
designed in the form of an annular rubber shear mount. A 
diagram of the damper is shown in figure 4.10 and is 
comprised of a brass casing filled with a soft prosthetic 
rubber of about 20 durometer hardness. A brass plunger was 
then cemented into the rubber in the center of the mount. An 
end of rod soldered to the casing was then screwed into the 
structure at point 19 while the opposite end of the plunger 
was screwed in at point 3. The design was based on the 
formulation for stiffness given by Macduff and Curreri [16] 
for rubber shear mounts. It was desired to minimize the 
amount of stiffness, considering physical constraints, and 
to accentuate the damping contribution. The longitudinal 
stiffness for annular rubber shear mounts 1S approximated 


ass 


Where bis the length, v5 and "| are the outside and inside 
radii respectively and G is the shear modulus for the 


rubber. 
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Figure 4.10 Annular Rubber Shear Mount Damper 
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A more accurate determination of the stiffness and damping 
was accomplished by testing the rod-rubber combination. This 
configuration was fixed with one end on a shaker table and a 
mass attached to the other end so as to approximate an 
oscillating support problem for a damped single DOF system. 
Accelerometers attached at the mass and the shaker ends 
allowed determination of the maximum ratio between the 
amplitude at the mass (X) and the amplitude at the shaker 


(Y) which for a linear viscous system forms the relation: 


where m is the mass on the end resonating at the input 
forcing frequency and c is the single unknown, the viscous 
damping constant. The stiffness of the rod is then 


calculated from: 


2 2 
k = MW ¢ 


The mass was selected so as to obtain a forcing frequency 
approximately midway in the frequency range of interest. The 
damping constant and stiffness were thus calculated to be 
0.13 lb-sec/in (23 N-sec/m) and 216 lb/in (3.78(10°*) N/m) 
respectively. 

The FRF, curve fit and coherence at point 1z for the 
damping addition to the structure is shown in figure 4.11. A 


relatively large amount of damping was expected as the time 
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Figure 4.11 FRF and Coherence for the Damping Modified 


Structure at Test Point 1z 
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Figure 4.12 Damping Modified Structure Synthesis Plot 
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response waveforms decayed in a fraction of the total time 
period. The high modal damping is again indicated by the 
wider peaks and obvious modal coupling as shown in the 
figure. nha the? cunve (hit iwas' relatively ‘good tatipcine 12, 
at most other points the fit was much worse. This difficulty 
BowciiplveetihcvDured. tovan limitation -ofmthus types orem. ele 
ESsPnotrknoOwnet other curvesfiitting algorithms sureer the 
Same limitations for high modal damping because of a lack of 
documentation in this area. This is perhaps due to the fact 
that most structures considered do not have such a high 
degree of damping. 

In modelling the change with SDM, it was again found 
that the additional mass and stiffnesses in the x and y 
directions played negligible roles in the dynamic behavior. 
It was assumed as well that damping in the x and y 
directions were also negligible for this type of mount. Thus 
the only modifications were those of damping and stiffness 
in the z direction between points 3 and 19. The compared 
results for the complex frequencies are shown in table 4.9. 
The predicted mode shapes are once again summarized in 
tables 4.10-4.14. The frequency test values from table 4.9 
for this modification are averaged values obtained at points 
where the curve fit appeared to conform closely to the 
experimental data. The table displays good agreement for the 
natural frequencies since they are most dependent on the 
stiffness modification, not the damping. The damping is much 


less consistent with a large discrepancy in the first mode. 
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In order to investigate the effects of possible variations 
in the experimental value of the damping constant, +20 % 
deviation in c was used to make predictions. These results, 
together with the original prediction and test results are 
shown in table 4.8. It can be seen that the damping in 
particular modes draws further away from the test values for 
both the higher and lower values of the damping coefficient. 
Therefore, it appears that the error in the damping 
coefficient cannot account for the effects discussed above. 

The synthesis plot for point 1z is shown in figure 
4,12. As can be seen, the differences between predicted and 
actual FRF amplitudes is much more than in the previous 
modifications. | 

The test structure and EMA data acquisition set-up are 


shown in plate 4.1. 
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MODE PARAMETER 


] FREQ (Hz) 
DAMP (2%) 
2 FREQ (Hz) 
DAMP (2%) 
3 FREQ (Hz) 
DAMP (%) 
4 FREQ (Hz) 
DAMP (3%) 
5 FREQ (Hz) 
DAMP (%) 
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Table 4.8 Effects on Prediction of a + 20% Variation in the 


Experimental Damping Constant 
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Table 4.12 Mode 3 Mode Shapes for the z Direction 
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Table 4.13 Mode 4 Mode Shapes for the z Direction 
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Table 4.14 Mode 5 Mode Shapes for the z Direction 
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Plate 4.1 EMA Data Acquisition Set-Up and Test Structure 
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5. CONCLUSIONS AND FURTHER STUDIES 

This study has shown how EMA and SDM can be combined to 
form a powerful tool in solving vibration problems. The 
prediction accuracy of modifications using SDM depends most 
highly on very good estimates of the original modal 
parameters from EMA since these form the modal model. The 
major inherent problems of noise and leakage contribute in 
reducing this accuracy so that the general case of complex 
modes is virtually impossible to accurately identify as 
such. As a result, a real mode model approximation was used 
where damping effects were taken into account through the 
assumption of proportional damping. Once the modal model had 
been constructed, accurate predictions using SDM were 
possible provided that the problems and limitations of the 
technique were realized. It was shown in chaper 3 how to 
reduce noise and leakage to improve the measurements. Noise 
greatly affects the measurements and thus predictions as 
displayed in the results obtained for the free-free beam. In 
using the modal model, it is important to realize that it is 
a truncated model with a frequency range set by the tests on 
fthesoriginalestructures yThus, ipredyvctions shar Youtside ‘the 
scope of the model could not be expected to be accurate as 
shown in the fifth mode prediction for the stiffness 
modification to the cantilever structure. Also found was 
that for many situations, especially for beamestructures, 
the rotary inertia effects become important as displayed in 
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these situations, rotational degrees-of-freedom must be 
estimated and added to the modal model. The final conclusion 
from this study was that damping effects at the best of 
times are difficult to predict since damping material 
properties typically vary a great deal, especially over 
frequency. In addition, the model used assumes proportional 
linear viscous damping. Therefore, damping modification 
predictions will ultimately be less accurate. 

Additional problems not included in this study but 
Still of considerable importance are those dealing with 
structural non-linearities and frequency resolution. When 
non-linearities are present in a structure, the excitation 
Spectrum level must be precisely controlled and thus the 
transient method becomes unsuitable. In this situation, the 
seme? classes of excitation, namely random and sinusoidal 
Signals are used. In the case of frequency resolution, the 
Zoom transform is employed to increase resolution by 
examination of discrete bands in the frequency range of 
interest. For the cantilever test structure used in this 
study, only frequencies up to 400 Hz were examined resulting 
in a very good frequency resolution of 1 Hz. However, in 
many situations, much higher frequencies are examined with a 
corresponding drop in frequency resolution and thus baseband 
resolution becomes inadequate. 

In determining a more accurate modal model, especially 
over high frequencies or when there is a high degree of 


modal coupling, Various improvements to the present 
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technique could be developed. These include the 
implementation of the Zoom transform and further 
investigation into improvement of the curve fitting 
algorithm. 

Further studies in the area of structural modifications 
might include modifications, especially damping 
modifications using a complex mode model like that described 
in O'Callahan and Avitabile [17] or using direct 
modification from FRF's as shown in reference [5]. 
Comparison of these techniques for various modifications and 
Ssubstructuring could yield further insight in this area. 

Extensions to uses of the modal model and SDM include 
forced response simulation and stress prediction. Both these 
techniques rely on the fact that the FRF completely 
describes a structures dynamics. Forced response simulation 
involves the determination of the response of the structure 
by multiplying the system FRF by a synthesized forcing 
function. Stress prediction combines the system FRF with 
Fourier transformed strain gauge operating time domain data 
to calculate a system 'pseudo' forcing function. A modified 
FRF found from using SDM can then be multiplied by the 
Noseudo “forcing fLumetion to obtail moditited system 
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